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*-£h \ The Next Linear Collider(NLC) will provide a excellent tool for probing the de- 

i ' 

Q_i' tailed nature of the top quark. By extending the recent analysis of Dokshitzer, Khoze 

<D ' 

and Sterling, we perform a preliminary examination of the influence of an anomalous 
■ chromomagnetic moment for the top, k, on the spectrum of gluon radiation associated 

with tt production. In particular, we analyze the sensitivity of future data to non-zero 
values of k and estimate the limits that can be placed on this parameter at the NLC 
with center of mass energies y/s = 500 and 1000 GeV. 
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Direct searches for the top quark at the Tevatron have led to a lower bound on the top 
mass of 131 GeV and, quite recently, to evidence that the top has been found with a mass 
near 175 GeV[|T], H. If verified, this will be a remarkable success for the Standard Model(SM) 
since this value of the mass lies close to the center of the range predicted by precision 
electroweak data||. Due to its large mass, the top itself has been proposed as a probe 
for new physics beyond the SM. Detailed analyses of top quark couplings to gauge bosons 
through its direct production and subsequent decay at both hadron[|J and e + e~ colliders|5| 
have been advocated in the literature for this very purpose. The present indirect constraints 
from low energy processes still allow for sizeable deviations from SM predictions @. In fact, 
the somewhat larger than expected cross section for ti production]?] obtained by the CDF 
Collaboration has already prompted several theoretical analyses || in which new dynamics 
involving the top quark have been discussed. Thus it is possible that the top may show us 
the first glimmer of new physics beyond the standard model. 

The possibility of gluon emission during heavy quark production in e + e~ collisions 
has been entertained for quite some time. In a recent paper, Dokshitzer, Khoze and 
Sterling(DKS)|§ have considered the spectrum of energetic gluon jets produced in asso- 
ciation with ti at the Next Linear Collider (NLC). In the present paper we will extend their 
analysis and consider the possibility that the top possesses a non-zero anomalous chromo- 
magnetic dipole moment, k, in its coupling to gluons. (As in the DKS study, we will ignore 
the effects of top decay in this analysis.) Such a scenario has recently been shown to lead to 



significant modifications in the characteristics of the ti production at the Tevatron |TI| . Here 
in this preliminary analysis we will show that the gluon energy distribution is a sensitive 
probe of k, as a non-zero value of this parameter leads to an enhancement in the number of 
high energy gluons produced along with ti at the NLC. Since the emission rate for gluons 
from a chromomagnetic dipole term in the Lagrangian scales approximately like K^fsjmt in 
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the amplitude, large values of k will lead to a breakdown in perturbation theory. We will 
also see that for a fixed top quark mass, increasing the NLC's center of mass energy will pro- 
vide an additional lever arm in obtaining sensitivity to non-zero values of k. In the absence 
of additional high energy gluon jets in comparison to ordinary QCD expectations we will 
demonstrate that reasonable limits on the value of k are obtainable at the NLC. One might 
expect that substantially improved limits from fits to the gluon energy spectrum itself may 
be obtainable. We show that vast improvements in the constraints on k from this approach 
are unlikely at the y/s = 500 GeV NLC due to a conspiracy between SM and ^-dependent 
contributions to the cross section. For the \/s = 1 TeV case, however, we find that fits to 
the gluon energy distribution yield greatly increased sensitivity to non-zero k. 

It is important to note that since the anomalous tig coupling does not occur at the 
primary production vertex at the NLC and appears only in a higher order process, the 
limits we expect on the value of k should be inferior to those obtained in the literature on 
anomalous electroweak top couplings ||. While this expectation is realized, the limits we 
obtain, particularly for the \/s = 1 TeV e + e" collider scenario, are reasonably strong. 

It is also important to remember, of course, that an 'anomalous' chromomagnetic 
moment for the top (or any quark) is induced at the one-loop level in conventional QCD and 
is of order 26 . In the context of the present paper, by 'anomalous' we mean a value over and 
above that given within the SM context, usually with a magnitude somewhat larger than 
that induced via conventional perturbative loop diagrams. 

To begin our analysis and in order to set our conventions, we note that the piece of 
the Lagrangian which governs the tig coupling with a non-zero value of k is given by : 
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where g s and T a are the usual SU(3) C coupling and generators, m t is the top quark mass, k 
is the out-going gluon momentum, and F 2 represents a general, k 2 — dependent, form factor. 
(A potential i<\(A; 2 )-type form-factor has already been set to unity.) For on-shell gluons, 
we define F 2 (k 2 — 0) = k following the standard notation. The rest of the notation below 
follows closely that of DKS. Let p±, p 2 , and k be the momenta of the t, i, and g in the final 
state such that q = p 1 + p 2 + k with q 2 = s. The kinematics of the e + e~ — > tig process then 
imply the usual defining relationships 

Zi = 2q- Pi/ s , 
z = 2q ■ k/s , 
2k-p 1 /s = I — Z2, 
2k-p 2 /s = l — z 1 , 

2m 2 

2 Pl -p 2 /s = 1-Z-—L, (2) 

where z 1 + z 2 + z = 2 is the statement of energy conservation. Following DKS, to leading 
order in a s we can factorize the weighted, angular integrated, double differential cross section 
for the process of interest into the separate contributions due to the vector and axial- vector 
couplings of the top quark to the s-channel exchange gauge bosons as 

d 2 W d 2 W v d 2 W a 

tr "~J J T-T a 3 1— , [o) 



dz\dz 2 v dz\dz 2 ° ' dz\dz 2 ' 

where F V;a are the 'weighting' factors telling us the fraction of events arising from the vector 
and axial-vector couplings of the top to 7 and Z. Note that we have scaled our result to 
the lowest order tt production cross section, i.e., W = a/a , where a = a(e + e~ — > tt). 
(Clearly, we must also have F v + F a = 1 in the above expression to conserve probability.) 
This factorization approach continues to remain valid even in the presence of a non-zero k. 



In order to proceed, only the quantities ^p^f and F v>a need to be computed. Since F v ^ a are 

insensitive to the existence of the anomalous chromomagnetic moment of the top quark, they 
are given solely by the kinematics and the electroweak couplings of the top. These factors 
can be read off directly from the general cross section expression for the production of heavy 
fermion pairs in e + e~ collisions, e.g., in ||ll|| : 



1/3(3 - p)A v 
v " ^A a + l(3(3-(3i)A v > 

F = U) 

^A a + 1/3(3- /3 2 )A' 1 1 



and /3 = y(l — 4mf/s), with s being the square of the center of mass energy. A v<a are 

directly determined by the vector and axial vector couplings of the electron and top quark 
as well as the gauge boson propagator function: 

A v = ( v i v i + "i'ljJ.ifir/JiPij , 

ij 

2 [(s - M?)(s - Mf) + (TM)i(TM)j] 
* " S [(a - Miy + (rM)f][( S - Mfy + {TM)]\ ' W 

The sum in the expression above is over the s-channel = 1) and Z(i,j = 2) gauge 

boson exchanges, including finite width effects, and for numerical purposes in our analysis 
the couplings are normalized to the running electromagnetic charge. In these numerical 
calculations, we will use the values of the various parameters as given in Ref . |J . We assume 
that the vector and axial vector couplings of the top quark are given by their conventional 
SM values with no alterations being present in the tij, Z vertices. 
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Given all of the above, the evaluation of the square of the matrix element is quite 
straightforward. Defining the overall normalization coefficients 

N v = ^(2m t V^)- 1 [^(3-^ 2 )]- 1 , 

N a = 2 ^(2mls 2 xjxl)-^r\ (6) 
where x« = 1 — Zi, we obtain the following expressions for the above distributions: 
d 2 W v 



N„ 



-8mj(a?i + x 2 ) - 4srrit[xj(l + 2x 2 ) + x 2 (l + 2x x )] + 2s 2 m 2 x 1 x 2 



dz\dz 2 

[(1 - xif + (1 - x 2 ) 2 + k(xi - x 2 ) 2 } + K 2 s 3 xlxl(l — x-i — x 2 ) 
d 2 W a r 

- — — — = N a 16m t (xi + x 2 ) 2 + 2sm\ [(k 2 + 2k + 2)x\x 2 {x 1 + x 2 ) 2 + 8xix 2 
dz\dz 2 L 

(xi + x 2 ) - 2(^1 + x\ + 6xix 2 )] + 2mjS 2 xix 2 [(l - £i) 2 + (1 - x 2 ) 2 + 
k(x\ + x\ — A) + k 2 xiX 2 (xi + x 2 — 3)] + k 2 s 3 x 2 x 2 (1 — Xi)(l — x 2 ) , (7) 

which, when combined with the other results above, gives the complete tig double differential 
cross section, normalized to that for ti production, including the contributions from finite k. 
For k = 0, we reproduce the standard results in the literature||. These expressions have no 
co-linear singularities, due to the finite top mass, but are still, overall, infrared singular in 
the limit of zero gluon momenta. It is already clear from Eqs. 6 and 7 that as s/ m 2 gets large 

the last terms in the expressions for will become dominant over the conventional QCD 

result. Note that these terms are infrared finite due to the fact that the chromomagnetic 
coupling in the Lagrangian is proportional to the gluon momenta. Semiquantitative bounds 
on k follow immediately upon integrating these two terms over phase space thus obtaining 
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the ratio of the ttg to tt cross sections in the large s/m 2 limit, i.e., 

(Jttg ^ a s k 2 s v 2 + 1.25a 2 

a a ~ vr 18m 2 v 2 + a 2 ' 1 j 

which leads for m t = 175 GeV and y/s = 500 GeV ( assuming that the large s/m 2 limit 
is crudely valid in this case) to ~ 5.0(k/3) 2 as the coefficient of — . Correcting for a color 
factor and a missing -, this limiting behavior is in agreement with that obtained by Grifols 



and Mendez[12] for an anomalous magnetic moment of the r in the decay Z — > rf^y. This 
crude estimate tells us that for perturbation theory to make sense for tig production at the 
NLC, the value of k must approximately satisfy \k\ < 3. 

There are several ways to see the effects of a non-zero k on the gluon jet energy 
spectrum. For the moment, let us follow DKS and calculate the average value of the scaled 
gluon energy, i.e., z ave , (where z = 2E g i u /^/~s) as a function of [3 2 for pure vector or pure 
axial vector couplings; these are shown in Figs, la and lb (assuming a s — 0.1 for purposes 
of direct comparison with DKS). Note that this quantity, z ave , is infrared safe for all values 
of (3. For k = 0, the results of DKS are immediately recovered but for non-zero k, significant 
deviations are observed which grow quite large with increasing f3 2 . It is, of course, just in 
this phase space regime where we expect the largest deviations from the conventional QCD 
results since large fl 2 implies large s/m 2 . Of course, even in standard QCD, the region close 
to (3 = 1 becomes non-perturbative (unless further cuts are applied) since it corresponds to 
the location in phase space where multiple soft co-linear gluon emission can occur. However, 
long before this non-perturbative range is reached we see that for non-zero values of k there 
is an upward shift in the average value of z, for both vector and axial vector couplings. This 
increased divergence in the average value of z for large f3 2 is only symptomatic of a more 
widespread phenomena, i.e., the the presence of a non-zero k not only increases the rate for 
the tig final state but also hardens the gluon jet energy distribution. We further note that 



the region of small (3 2 is also nonperturbative as this corresponds to the threshold regime 
where resummation techniques need to be applied. The possible signatures for a non-zero k 
from detailed threshold region studies is beyond the scope of the present paper. 

Perhaps a better probe of non-zero k is the gluon jet energy spectrum itself. This we 
show in Fig. 2a assuming mt = 175 GeV and an NLC center of mass energy of \/s = 500 
GeV. From the kinematics, the maximum gluon energy for these input parameters is z max = 
1 — 4mf/s = 0.51 or 127.5 GeV. Note that in the Figure we have again assumed a s = 0.1. 
To get the 'best' result to compare with experiment (ignoring higher order corrections) we 
need to obtain the correct value of a s at this energy and simply rescale all of the curves in 
Fig. 2a by an overall factor. One possible approach, which we will follow below, is to use 
the method of Brodsky, LePage, and Mackenzie(BLM) ||13|| in order to set the scale, using 
the value of a s (M 2 )0 as input, and to make use of the three-loop renormalization group 
equations. Without doing a complete calculation, it is possible to estimate the correct BLM 
scale(<5*) at which to evaluate a s by simple phase space considerations WM. One finds 



Q* ~ 0.435( v /s — 2m t )/3, not too different than the Q* for Z decay studies, so that we arrive 
at a s (Q*) = 0.121 as our estimate (using the value of a s at the Z scale extracted from -R/i|| 
as input). This implies that all the curves in Fig. 2 should be scaled upward by about 21%. 
For other preferred values of a s , an appropriate rescaling should be performed. 

Of course, the overall scale is not the interesting part of these figures as far as k sen- 
sitivity is concerned since k non-zero results in a distortion in the gluon spectrum, especially 
at larger z values. Generally one sees that the effect of k is to flatten the spectrum so that 
there is an excess of gluon jets with high energies. For k = —1, however, we see a sharper 
fall off in the spectrum than in standard QCD; this is a result of a destructive interference 
between the ordinary and K-dependent amplitudes which takes place for the specific values of 
m t and y/s we have used as input. This implies that for a range of negative k, the SM result 
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and the one where k is non-zero will be very difficult to distinguish. Since the deviation 
due to k is at larger z values and the spectrum diverges as z — > 0, we can apply a cut on 
the minimum gluon energy, z cut , which we use to define our event sample and integrate the 
spectrum for gluon energies above that value. This results in Fig. 2b, which should also be 
scaled upwards by 21% if the BLM approach is used. As expected, the curves for non-zero k 
are generally higher than the standard QCD result. Taking the BLM value for a s and a value 
of z cut = 0.2, i.e., only events with gluon jets having energies in excess of 50 GeV, we show in 
Fig. 2c the k dependence of the resulting integrated cross section. Assuming an integrated 
luminosity of SOfb^ 1 , this corresponds to a sample of 375(7/0.02) tig events before further 
cuts are applied or 457 events in the standard QCD case with k — 0. Of course to identify 
top pair production, we will demand at least one high-p t lepton in the event (B = 0.44) and, 
perhaps, an additional &-tag(with an assumed efficiency of e = 0.8) to remove backgrounds. 

If the SM result is realized, we can use the estimates of the tig event rate above to place 
bounds on the value of k. (We might expect better limits would most likely be obtainable 
by a direct measurement of the gluon energy spectrum instead of a simple rate estimate; we 
will return to this possibility below.) Allowing for a 2% error on the determination of a s 
in the NLC era, a 5% systematic error from higher-order QCD uncertainties, and the rates 
above to calculate statistical errors, the value of I would be determined to be 2.44 ± 0.23%, 
which at 95% CL would restrict k to lie in the range —2.1 < k < 0.6. Varying the cuts leads 
to numerically similar results there being a relative trade off between increased(decreased) 
statistics and decreased(increased) sensitivity. The reason for the poor limit is clear; as we 
saw above a range of negative k exists for which the resulting energy distribution is almost 
identical to the SM one. This is a result of a conspiracy between the various contributions 
to the differential cross section and the particular values of m t and y/s we are examining. 

This situation is not much alleviated by an actual fit to the gluon energy distribution 
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itself. For z > 0.15, we generate 'data' assuming k = by Monte Carlo taking a bin size 
of Az = 0.05. The phase space region of interest is then covered by a total of seven z bins, 
the last one covering the range 0.45 < z < (3 2 = 0.51. We then fit the ^-dependent gluon 
spectrum to the data and perform a x 2 analysis. The resulting x 2 is a quartic function of 
k and is shown explicitly in Fig. 2d where one sees that two essentially degenerate minima 
exists. This is due to the conspiracy discussed above and, as a result, only rather poor 
bounds on k are obtainable. Explicitly, from this procedure we obtain the allowed range 
— 1.98 < k < 0.44 at 95% CL, which is only a slight improvement in the limit obtained 
above from simple counting. Similar limits are obtained if different bin sizes are chosen for 
the Monte Carlo study but we have not tried to optimize this choice in our analysis. In Fig. 
2e, we compare the data generated by Monte Carlo with the energy spectra predicted for 
the k values corresponding to the two approximately degenerate x 2 minima, i.e., k = —1.69 
and 0.12. As can be seen the fit is quite good in both cases. 

What happens when we go to larger values of y/s where significantly greater sensitiv- 
ities to k are expected? In Figs. 3a-d, we examine the case of tig production with non-zero 
k for an NLC with y/s = 1 TeV. (In these figures, the BLM value of a s = 0.100 has been 
assumed following the above procedure so that no overall rescaling is necessary in this case.) 
Fig. 3a clearly shows that at large values of z, the gluon jet energy distribution is even 
more enhanced for fixed values of k than for the yfs = 500 GeV case; this is exactly what 
we should have expected. Note the approximate symmetry of the curves under the inter- 
change k — > —k; this occurs naturally in the large s/mf limit as seen above. Applying the 
z C ut approach as before yields Fig. 3b where the k — > — k symmetry is even more obvious. 
The value of / is so large in the \k\ — 3 case for values of z cut < 0.4 that we should most 
likely not trust our lowest order perturbative result for this range of parameters. Taking 
z C ut = 0.4, which corresponds to a gluon jet of energy 200 GeV, we show in Fig. 3c the 
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explicit re dependence of our cross section ratio. If we assume an integrated luminosity of 
200/6" 1 and make the same assumptions as in the 500 GeV case, the realization of the SM 
result can again be used to place significant constraints on re. In this case we would obtain 
I = 5.24 ± 0.36%, which at 95% CL would restrict re to lie in the range -1.0 < re < 0.25, 
which is about a factor of two better than that obtained for the 500 GeV NLC. We see again 
that because small, negative values of re and the SM case lead to very similar gluon energy 
distributions, we do not obtain a very narrow allowed range for re. However, if we try the 
Monte Carlo approach as we did above and fit the entire spectrum for z > 0.4 in nine z bins 
(the last bin covering the range 0.8 < z < j3 2 = 0.8775) we obtain the x 2 plot shown in Fig. 
3d. Unlike the t/s = 500 GeV case, the second local x 2 minima is no longer degenerate so 
that we can now obtain a substantially improved bound on re: —0.12 < re < 0.21. Again, 
we have not made any attempt to optimize the bin size in this Monte Carlo study. In Fig. 
3e, we compare the Monte Carlo generated data with the predicted gluon spectrum for the 
choice of re = 0.06, corresponding to the x 2 minimum. As can be seen, the fit is quite good. 

In this paper we have analyzed the influence of an anomalous chromomagnetic dipole 
moment for the top quark, re, on the cross sections and associated gluon jet energy distri- 
butions for tig events produced at both 500 and 1000 GeV e + e~ linear colliders assuming 
a top quark mass of 175 GeV. Making a cut on the gluon jet energy of 50(200) GeV and 
demanding at least one b-tag as well as one high-p t lepton tag as a top signal, an integrated 
luminosity of 30(200) fb' 1 leads to a bound on re of -2.1 < re < 0.6(-1.0 < re < 0.25) 
at the NLC with y/s = 500(1000) GeV, assuming no excess tig events are observed. One 
might have expected that a complete fit of the gluon energy spectrum to the re-dependent 
distribution would generally lead to substantial improvements in these limits. We found, 
however, that this was not the case at a \/s = 500 GeV NLC due to a conspiracy between 
the SM and re-dependent terms in the cross section and that only slight improvements were 
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obtainable: —1.98 < k < 0.44. For the case of a y/s = 1 TeV machine, this conspiracy no 
longer took place and the power of fitting to the gluon energy spectrum was realized yielding 
a vastly improved bound of —0.12 < k < 0.21. 

The results of this analysis are only preliminary. In a more complete Monte Carlo 
study the effects of the top decay possible gluon emission from the final state bottom quarks, 
and detector resolution and efficiencies need to be included. However, by demanding a very 
high energy additional gluon jet with the rest of the event reconstructing to ti, the results of 
such an analysis should closely mimic those we have obtained above. We have seen that it is 
quite likely that the NLC will be able to place reasonably strong constraints on the existence 
of a top quark anomalous chromomagnetic moment. Of course, the more exciting possibility 
of observing a hardening of the gluon jet energy spectrum would be a spectacular signature 
for an anomalous magnetic moment for the top or other new physics beyond the standard 
model. 
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Figure Captions 

Figure 1. Average value of the scaled gluon jet energy, z, as a function of f3 2 for both purely 
vector(a) or axial vector(b) heavy quark couplings. The upper(lower) dotted, dashed, 
and dot-dashed curves correspond to k values of 3(-3), 2(-2), and 1(-1) respectively 
while the solid curve is conventional QCD with k — 0. a s = 0.10 has been assumed. 

Figure 2. (a)Gluon jet energy spectrum assuming a s = 0.10 for m t = 175 GeV at a y/s = 500 
GeV NLC. (b) Integrated gluon energy spectrum for the same input parameters as 
in (a) as a function of the minimum gluon energy, z cut . In both (a) and (b), the 
labelling of the various curves is as in Fig. 1. (c) Same as (b) but as a function of k 
assuming z cut = 0.2 and a s = 0.121 as suggested by the BLM approach, (d) x 2 fit to 
the gluon energy spectrum for the value of k as described in the text. a s = 0.121 is 
again assumed, (e) Best fit gluon ^-dependent spectra through the points generated 
by Monte Carlo. The dashed(dash-dotted) curve corresponds to to k — —1.69(0.12). 

Figure 3. Same as Fig. 2, but for an NLC with y/s = 1 TeV. In (c) the value of z cut = 0.4 
is assumed along with a s = 0.100. (e) Best fit gluon spectrum through the points 
generated by the Monte Carlo analysis corresponding to k — 0.06. 
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